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Abstract 



We will give a review of the computation of exact next-to-leading order corrections to 
heavy quark production in deep inelastic lepton-hadron scattering and discuss the progress 
made in this held over the past ten years. In this approach, hereafter called EXACT, 
where the heavy quark mass is taken to be of the same order of magnitude as the other 
large scales in the process, one can apply perturbation theory in all orders of the strong 
coupling constant a s . The results are compared with another approach, called the variable 
flavor number scheme (VFNS), where the heavy quark is also treated as a massless quark. 
It turns out that the differences between the two approaches are very small provided both 
,^ ' of them are carried out up to next-to-next-to-leading order. 
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£2 1 Introduction 



In the last ten years one has made much progress on the theoretical and experimental level in 
the study of heavy flavor production in deep inelastic lepton-hadron scattering. Computations 
of the cross section in the Born approximation, where the heavy quark is treated as a massive 
particle, were finished by the end of the seventies [0 whereas the first experimental results came 
from the EMC-collaboration |J in the early eighties. In the nineties the 0(a s ) corrections to the 
Born process were computed in ||]. Moreover other methods to study heavy flavor production 
were advocated like the intrinsic quark approach j|] and the variable flavor number scheme || 
where the heavy quark is also treated as a massless particle. From the experimental side the 
electron (positron)-proton collider HERA has given us a wealth of information about charm 
quark production (for recent experimental results see ||, [[/J, ||, 0). We will report about this 
progress below. 

The semi-inclusive process describing heavy quark production in deep inelastic electron- 
proton scattering is given by (see Fig. [Xj) 

e-(fci) + P(p) -> e -(fe) + Q( Pl ) +' X' , (1) 
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Figure 1: Kinematics of heavy quark Q-production in deep inelastic electron-proton scattering 



where V represents the intermediate vector boson Z or 7 carrying the momentum q and Q 
denotes the heavy (anti-) quark. The scaling variables are given by 
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In the ongoing experiments carried out at HERA |J, @, ||, [|] and in the fixed target experi- 
ments carried out in the past 0, Q 2 <C M§ so that the intermediate Z-boson can be neglected 
and the reaction is dominated by photon exchange. In this case the inclusive cross section 
simplifies considerably and when the incoming particles are unpolarised it can be written as 
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where F^q (k = 2, L) denote the heavy quark contributions to the structure functions. Notice 
that an analogous formula exists for the semi-inclusive cross section 



d A a 



dx dy dp TtQ dr] Q 
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where Pt,q and i]q denote the transverse momentum and rapidity of the heavy quark Q respec- 
tively. Both of them are considered in the center of mass of the photon-proton system. For 
neutral current reactions one has proposed two different production mechanisms for heavy quark 
production. They are distinguished as follows 
I Intrinsic Heavy Quark Production 

Here one assumes that, besides light quarks u, d, s and gluons g, the wave function of the 
proton also consists of heavy quarks like c, b, t ||. In the context of the QCD improved 
parton model this means that the production mechanism is described as indicated in Fig. 
H]. In this picture the heavy quark emerges directly from the proton and interacts with 
the virtual photon 7*. The consequence is that it is described by a heavy flavor density 
/g(z,/i 2 ) with pg = zp and fi denotes the factorization scale. For this mechanism the 
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Figure 2: (a) 7* + Q -> Q (#£), (6) 7* + Q - Q + 9 (4%). 

heavy quark structure function has the following representation 

/■l dz 

F k ,Q(x,Q 2 ,m 2 ) = e 2 Q / — f Q (x/z, fi 2 ) H k , Q (z,Q 2 ,m 2 , fi 2 ) 

J x Z 

= e 2 Q f Q ^ 2 )®H KQ {Q 2 ,m 2 ^ 2 ). (5) 

Here eg denotes the charge of the heavy quark and ® stands for the convolution 

/®^(x)= / -f(x/z)g{z). (6) 

Further the heavy quark coefficient function i/^Q can be expanded in the strong coupling 
constant a s (/i 2 ) as follows 

00 

tfw =£«:#$, ^th 0. = /. (?) 

n=0 47r 

Some of the contributions to -f/^Q are given by the diagrams in Fig. |]a,b. 

II Extrinsic Heavy Quark Production 

In this case the proton wave function does not contain the heavy quark components. In 
lowest order of perturbation theory the heavy quark and heavy anti-quark appear in pairs 
and are produced via photon-gluon fusion |lj as presented in Fig. ||. Here the gluon 
emerges from the proton in a similar way as the heavy quark in Fig. 0. In this approach 
the heavy quark structure function reads in lowest order 

F k , Q (x,Q 2 ,m 2 ) = a s e 2 Q f g ^ 2 ) ® tf«(Q 2 ,m 2 ,/i 2 ) . (8) 

The main difference between the two production mechanisms can be attributed to the fact that 
for extrinsic heavy quark production two heavy particles are produced in the final state instead 
of one as in the case of the intrinsic heavy quark approach. This reveals itself in the transverse 
momentum pr-distribution where for mechanism II the quark and anti-quark appear back to 
back. The experiments carried out HERA J9J confirm the pr-spectrum predicted by the latter 
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Figure 3: Feynman diagrams for the lowest-order photon-gluon fusion process contributing to 
the coefficient functions H\. g . 



mechanism. However in the past the EMC-collaboration []2[ found a discrepancy at large x 
(x = 0.237) between the structure function F 2 ^q, predicted by mechanism II, and the charm 
quark data. This difference can be explained by also invoking mechanism I || (see also pi ) . 
Nevertheless because of the success of extrinsic heavy quark production revealed by the HERA 
charm quark data we will continue with this approach in our calculations below. 



2 Heavy Quark Structure Functions up to 0(a 



In the extrinsic heavy quark approach there are two different production mechanisms P"T| . 
T2| which are of importance for the derivation of the variable flavor number scheme (VFNS) 
treated in the next section. In the case of electro-production the virtual photon either interacts 
with the heavy quark appearing in the final state or it is attached to the light (anti-) quark in 
the initial state. This distinction is revealed by the Feynman graphs. Some of them will be 
shown as illustration below. 

A The photon interacts with the heavy quark 

Here the lowest order (LO) process is given by gluon fusion as presented in Fig. |3|. It is 
given by the partonic reaction 

1*+9^Q + Q, (9) 

which is calculated in 0. In next-to-leading order (NLO) one has to compute the virtual 
corrections to this process. Some of the Feynman graphs are shown in Fig. |j. Besides the 
virtual corrections one also has to calculate all processes which contain an extra particle 
in the final state. Adding a gluon to reaction @ we observe gluon bremsstrahlung given 
by 

l* + 9^Q + Q + 9- (10) 

Some of the Feynman graphs corresponding to the process above are presented in Fig. [5]. 
Besides this reaction we have another one which appears for the first time if the computa- 
tions are extended to NLO. It is represented by the Bethe-Heitler process 

i* + q(q)^Q + Q + q(q), (n) 




Figure 4: Virtual gluon corrections to the process 7* +g — > Q + Q contributing to the coefficient 
functions iff) 
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Figure 5: The bremsstrahlung process 7* +g — > Q + Q+g contributing to the coefficient functions 

h9 



where the Feynman diagrams are shown in Fig. §. The contributions to the heavy quark 
structure functions due to the reactions above can be written as 
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(12) 



Here eg denotes the charge of Q indicating that the photon couples to the heavy quark and 
Hk,i (i = q,g), characteristic of these type of processes, represent the heavy quark coefficient 
functions which emerge from the calculation of the Feynman graphs after renormalization 
and mass factorization is carried out. The former procedure leads to a dependence on a 
renormalization scale \x of the coupling constant a s , the coefficient function H^i and the 
parton densities /$. Moreover the latter two also depend on a factorization scale which 




Figure 6: The Bethe-Heitler process 7* + <3 — > Q + Q + q contributing to the coefficient functions 

rr(2) 




vwvwv- 



-TODDK^ 



Figure 7: The Compton process 7* + Q — > Q + Q + q contributing to the coefficient functions 

r(2) 



for convenience is set equal to the parameter \x defined above. The factorization scale, 
which can be attributed to mass factorization, appears in all coefficient functions except 
for the lowest order one i.e. H ki . In the structure function of Eq. (fTj ) there appear two 



different types of flavor singlet parton densities i.e. the gluon density f g and the quark 
singlet density which in a three flavor number scheme reads as 

/ 9 S (3,/i 2 ) = 

/«(3, /i 2 ) + /«(3, fi 2 ) + f d (3, /i 2 ) + / d -(3, /i 2 ) + / a (3, /i 2 ) + f s (3, tf) . (13) 

Hence Eq. (|T2p represents the singlet part of the heavy quark structure function. 
B The photon interacts with the light quark. 

The second production mechanism is represented by the Feynman diagrams where the 
photon couples to the light quark or the anti-quark. This happens for the first time in 
NLO where one observes the Compton process 

Y + q(q)^Q + Q + q(q), (14) 

which is depicted in Fig. [F| The coefficient functions corresponding to this type of process 
are denoted by L ki (i = q, g) and the contribution to the heavy quark structure functions 
is characterized by the expression 

F fe , Q (Q 2 ,m 2 ) = £ ^oVlf/A^ + ^^g^m 8 ), (15) 

i=u,d,s 

where e« denotes the charge of the light quark represented by i = u, d, s in a three fla- 
vor number scheme. Since this process appears for the first time in second order mass 

(2) 

factorization is not needed which explains the independence of L k on the parameter \i. 

The computation of the second order contributions to the heavy quark coefficient functions 
H k ,i, L kt i has been carried out in 0. While calculating the Feynman graphs in Figs. £§- 1| one 
encounters several type of singularities which have to be regularized and subsequently to be 
subtracted off before one obtains a finite result. The singularities are of the following nature 
i.e. infrared (IR), ultraviolet (UV) and collinear (C). Sometimes the latter are also called mass 
singularities. The IR divergences cancel between the virtual and the bremsstrahlung corrections 
to reaction (|9|). The UV divergences, regularized by n-dimensional regularization, are removed 
by mass and coupling constant renormalization. For the mass renormalization we choose the 



on-shell scheme. In this case the UV divergence will be removed by replacing the bare mass m 
by the renormalized mass m via 
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where the UV pole term is indicated by 1/e with e = n — 4. If we choose for example process ( |Tol) 
together with the virtual corrections to the Born reaction (||) (see Figs. |],[|) the unrenormalized 
coefficient function takes the form 
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where H k ^\ fJ-=m is finite and d s denotes the bare coupling constant. Further we have also regular- 
ized the collinear divergences by n-dimensional regularization. In order to distinguish between 
ultraviolet and collinear divergences we have indicated them by 1/euv and 1/ec respectively. 
The residues of the collinear divergences are represented by the so called splitting functions de- 
noted by Pij (i,j = q,g). The origin of the collinear divergences is explained by the first diagram 
in Fig. H The propagator carrying the momentum k\ — k 2 behaves as 



{k\ — k 2 ) 2 2cg>! lu 2 (1 — cos#) ' 

which diverges for 9 — > 0. This propagator only shows up if three massless particles are coupled 
to each other like three gluons or when a gluon is attached to a quark line provided the quark 
is massless. If the gluon is attached to a heavy quark the mass of the latter, which is unequal 
to zero, prevents that the denominator in Eq. (0) vanishes when 6^0. The UV-divergence in 
Eq. (|i~7l) is removed by coupling constant renormalization which is achieved by adding a s H\, 



to Eq. 
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T7\) and replacing the bare coupling constant a s by the renormalized one represented by 
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(19) 



Finally one has to remove the collinear divergences. This is achieved by mass factorization. It 
proceeds in a similar way as multiplicative renormalization so that one can write 
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where T^ represents the transition function which removes all collinear divergences from the 
bare heavy quark coefficient functions. Further we have introduced the notion of bare parton 
density /j so that the heavy quark structure function can be written as 
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F k ,Q(Q , m ) = e I Q f l ® H Ki — , Q\ m' 



Substitution of H ki (see Eq. (P0|) ) into the expression above we can derive 

F fc , Q (Q 2 ,m 2 ) = <? Q fj{n 2 )®H h j (Q 2 ,m 2 ^ 2 " 



(22) 



where fj(fi 2 ) is the renormalized parton density defined by 



/ 3 V) = r i< (— ,^)®/i. 



(23) 



In the case of the example presented in Eq. (pl|) the collinear divergence is removed by adding 
H^' to H^g and choosing the following transition function 
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where N denotes the number of colors. 






C? 



c 



b 



10 
9 

8 
7 
6 
5 
4 
3 
2 
1 




i 1 1 1 r 

Eur. Phys. J. C12 (2000) 35: ZEUS 96-97 I — h 

Osaka Abs. 855: ZEUS 95-00 ! — h 

Osaka Abs. 984: HI 96-97 :- - X- 

EXACT 

BMSN 

CSN -■-■ 



Jr.-*---; 



. hill 






^fi 



^*a 



P*3' v 






0.5 



1.5 



logio Q 2 



Figure 8: The combined Osaka and published Zeus-data for da / dlog^Q 2 in nb for deep inelastic 
production of Z)* ± -mesons. The dashed line is the exact NLO result from the program HVQDIS 
which follows from F 2 E ^ CACT . The dotted line (BMSN-scheme) and dashed-dotted line (CSN- 
scheme) is based on F^f NS . 

A comparison of the next-to-leading order (NLO) heavy quark structure function F 2tC with 
the data for charm quark production measured at HERA reveals a fairly good agreement between 
theory and experiment. The data cover the range 1 < Q 2 < 1350 GeV 2 and 5 x 10~ 5 < x < 
5.6 x 10~ 2 . In [[U| one has made a comparison with the data obtained by the ZEUS-collaboration 
0,0. From the cross section in Eq. (0) one can derive the integrated quantities 
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(25) 



Notice that the quantities above represent .D*-meson production rather than charm quark pro- 
duction. The meson appears as a fragmentation product of the quark and the cross sections in 
Eq. ( p5|) are obtained by convoluting Eq. (|2|) with fragmentation functions. Furthermore one 
has to impose experimental cuts on the kinematics which are indicated by max and min. The 
results are presented in Figs. §, |9| which originate from [[U| where one can also find the maximal 
and minimal values for x and Q 2 . The figures show that there is a good agreement between 
the exact NLO result (called EXACT in the figure) and the data except for x ~ 1CT 3 where 



there is a small discrepancy. Furthermore in [jy|] one has also made a comparison between the 
program HQVDIS [|HJ based on the NLO computations above and the experimental differential 
distributions for D* production where the following cross sections are studied. 
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where W (W 2 = (p + q) 2 ) is the center of mass energy of the photon-proton system. Also the 
differential distributions agree with the NLO predictions except for the rapidity ^-distribution. 
Here it appears that for rj D > the experimental cross section is larger than the one computed 
in NLO (see also IB). 
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Figure 9: The combined Osaka and published Zeus-data for da / 'dlogiQX in nb for deep inelastic 
production of D* ± -mesons. The dashed line is the exact NLO result from the program HVQDIS 
which follows from F 2 E ^ ACT . The dotted line (BMSN-scheme) and dashed-dotted line (CSN- 
scheme) is based on F^, J NS - 



Summarizing our findings for the exact NLO calculations we conclude 
1. There exist a fairly good agreement between the data and the NLO calculations. 



2. The theoretical curves are insensitive to the choice of the renormalization/factorization 
scale \i occurring in the parton densities and the coefficient functions 



3. The theoretical curves are sensitive to the choice of the charm mass m c which is in the 
range 1.3 < m c < 1.7 GeV/c 2 



~2 



4. Processes with a gluon in the initial state (Eqs. (P),(|TUD) constitute the bulk of the con- 
tribution to the heavy quark structure function at x < 10 -2 . Hence they are an excellent 
probe to measure the gluon density f g (z,n 2 ). 

3 Asymptotic Heavy Quark Structure Functions 

In this section we want to study the heavy quark coefficient functions Hk.i, Lk^ in the asymp- 
totic region Q 2 3> m 2 . In this region they have the following form 

ff#™(z, Q 2 , m 2 , f) ~ £ a\ £ a nm {z) In" (£) In™ g) . (27) 

1=1 n+m<l \ ln / \ ln I 

A similar expression exists for L^f YMP . An example is the asymptotic expression for the Born 
reaction in Eq. (S) (Fig. |) which can be written as 



H^™ p (z,Q 2 ,m 2 )=a s 



l (Q 2 \ 

-P qg {z) In I — j + a Qg {z) + c 2 , g (z) 



(28) 



The origin of this asymptotic behaviour can be attributed to the property that in the limit 
m — > the heavy quark coefficient functions become collinear divergent which is revealed by the 
logarithmic singularities InQ 2 /m 2 and ln// 2 /m 2 . The reason why this behaviour is of interest 
can be summarized as follows 

1. The results obtained for the exact coefficient functions in the previous section are semi- 
analytic. In one has obtained exact results for the virtual corrections to process (§) 
(Fig. |) but for the reactions with three particles in the final state like Eqs. JlOf) , (PT) , 
( PI) a full analytical expression could only be presented for the Compton process in Eq. 
(0) ( see [['3D- m the other cases only the integration over the angles could be carried 
out but the integration over the final state energies are so tedious that they have to be 
done in a numerical way. However the latter integration becomes more amenable when 
m 2 <C Q 2 so that terms of the order m 2 /Q 2 can be neglected. Therefore an analytical 
result for the asymptotic heavy quark coefficient functions provides us with a check of the 
exact expressions computed for arbitrary m. 

2. The asymptotic coefficient functions play an important role in the derivation of the variable 
flavor number scheme ||, [17|] discussed at the end of this section. This scheme is only 



useful if the following questions are answered. They are: 

a. Do the logarithmic terms of the type ln n Q 2 /m 2 , occurring in the coefficient functions, 
really dominate the heavy quark structure functions or the heavy quark cross sections? 

b. Are the logarithmic terms ln n Q 2 /m 2 so large that they bedevil the perturbation series 
so that they have to be re-summed? 

There are two ways to compute the asymptotic heavy quark coefficient functions. 
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1. One can follow the procedure for the derivation of the exact coefficient functions in || but 
since the mass m can be neglected one can now carry out the additional integrations over 
the energies of the final state particles in an analytical way. 

2. One can use the operator product expansion (OPE) techniques which however are only 
applicable for inclusive quantities like the structure functions F kt Q(x,Q 2 ,m 2 ). 



In |TT[ one has adopted the latter approach which will be outlined below. In the derivation of 
the cross section in Eq. (|3|) one encounters the hadronic tensor leading to the definition of the 
structure functions. It is defined by 
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(29) 



In the limit Q 2 ^> M 2 the current-current correlation function appearing in the integrand of Eq. 
( PP| ) is dominated by the light cone. Hence one can make an operator product expansion near 
the light cone and write 

J{*),m] ~ EE^Wv-^of-^to^ 2 ), (so) 



where for convenience we have dropped the Lorentz indices of the currents. Here c 4 denote the 
coefficient functions which are distributions and Off 1 '"* 1 ™ are local operators. Both are renormal- 
ized which is indicated by the renormalization scale p,. When dropping all terms of the order 
M 2 /Q 2 we can limit ourselves to leading twist operators. In QCD they are given by 
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(31) 



The symbol S indicates that one has to symmetrize over all Lorentz indices. The covariant 
derivative is given by D^ = d^ — igT a A^ and A r is the generator of the flavor group SU(nf)p. 
Insertion of the OPE (Eq. (0)) into the hadronic tensor of Eq. ( |29"D leads to the following result 
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where the operator matrix element and the coefficient function are defined by 

4 M (^) = (P( P )\Or(0,v 2 )\P(p)) , Ct ] {^Pj = j*z<*'<t\*f) . (33) 

The OPE techniques can also be applied when the the proton state \P(p)) in Eq. (|29| ) is replaced 
by a light quark state \q{p)) or a gluon state \g(p)). However when the proton is replaced 
by massless quarks and gluons the external momentum satisfies the relation p 2 = so that 
the partonic structure functions and the partonic operator matrix elements become collinearly 
divergent. One can show (see [I7|) that instead of Eq. ( ft2| ) one obtains more complicate 
expressions which are given by 

f-NS I Q_ J_ | _,_ fASYMP I Q_ ^_ J_ | JNS ( m J_ 



^U 2 '^ M U'V'ecJ w U 2 ^ 




r-S | V ' I _i_ f ASYMP I ^ m ' 1 i rVASYMP (V m 

,Afc ' i I -,2 ' c„ + k * Un2' ,,2 ' 7" +J "fe,i Ur.2' „2 



s I 2L_ __ I _|_ fASYMP / ^t_ 'JJ_ J_ \ | fVASYMP 

/i 2 'e c / *'* lm 2 '/i 2 'W k/l \m* ^ e c 



iJ = Q,9- ( 34 ) 

Here !Fk,i represent the partonic structure functions which are given by Feynman graphs contain- 
ing massless particles only and therefore contain collinear singularities indicated by ec- These 
singularities also appear in the asymptotic heavy quark coefficient functions H^f YMF , L^f YMP 
which are determined in the asymptotic regime Q 2 ^> m 2 before mass factorization is carried 
out. These collinear divergences can be traced back to the massless quarks and gluons appearing 
in Figs. 0-0. Finally Aji represent the operator matrix elements on the partonic level and are 
defined by 

/™2 1 \ 

(i\Oj(jjL 2 ,Q)\i), i = q,9, J = q,9,Q, P 2 = 0. (35) 

The operator matrix elements which depend on the heavy flavor mass consist of two classes. The 
first class is given by heavy quark operators sandwiched between gluon or light quark states. 
The second class contains light quark or gluon operators which contain a heavy flavor loop. 
Examples of the first class together with the corresponding process are given in Figs. [Tl] and [TTJ. 
An example of the second class is shown in [12]. The light partonic coefficient functions defined 
by Ckj are derived from the light partonic structure functions via mass factorization 

(36) 






where r^ denote the transition functions which are discussed below Eq. (2C). The quantities 



J-'k ,i and Ck,i have been calculated up to second order in a s in ]nj. The operator matrix elements 
Aij are also known up to second order and the calculation of them is presented in [ 11 1 , [ 17] . Like 
in the case of the coefficient functions n-dimensional regularization has been used to regularize 
the ultraviolet and collinear divergences and one has chosen the same renormalization conditions 
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Figure 10: The operator matrix element A)v! and the corresponding Feynman graph for the 
process 7* + g -> Q + Q (HP). 





(2) 
Figure 11: The operator matrix element Aq and the corresponding Feynman graph for the 

process Y + g^Q + Q + g (H^ g ). 



13 



/ \ 

/ \ 

/ \ 



/ Q \ 

/ Q \ 



7 



g* 



/ 



Q 

Q 



7 



h 



jms\ 



\ 



Figure 12: The operator matrix element A^'^ and the corresponding Feynman graph for the 
process Y + Q^Q + Q + Q ( L ^l )■ 



for the mass and the strong coupling constant. Hence from Eq. (|34|) one infers the asymptotic 
heavy quark coefficient functions H^f YMP , L^f YMP up to the same order (see |T7[] ). Finally one 
can remove the remaining collinear divergences via the same mass factorization as is done for 
the exact heavy quark coefficient functions in Eq. (PUD- After this outline of the calculation 
of the asymptotic heavy quark coefficient function one might ask the question whether it is 
not easier to compute them in a more direct way. The main problem of radiative corrections 
is the computation of the phase space integrals in particular if one has massive particles in 
the final state even if one takes m? -C Q 2 . Since this work was already done for the light 
partonic structure functions Tk,% in [[H| it was not needed to repeat this procedure anymore. 
On the contrary it is much easier to compute two-loop operator matrix elements because of the 
zero momentum flowing into the operator vertex indicated by the symbol 
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in Figs. 

The difference between the exact and asymptotic coefficient functions can be attributed to the 
power corrections of the type (m 2 /Q 2 ) 1 ,with I > 1, which occur in the former but are absent 
in the latter. The asymptotic coefficient functions only contain the logarithms In" 1 Q 2 /m 2 and 
ln n n 2 /m 2 and terms which survive in the limit Q 2 — > oo. We will now first answer the question 
whether these logarithmic terms dominate the heavy quark structure function Ff,^Q(x,Q 2 



For that purpose we have studied the structure functions for charm production i.e. 
TTfl , fT9|j . Here one has computed the ratio 



Q 



m ). 
■ c in 



R 



k.c 



pASYMP 
pEXACT 



(37) 



where F^ and F^ YMP are represented in the three flavor number scheme 
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Figure 13: i?2,c i n NLO (Eq. fl37|) ) plotted as function of Q 2 at fixed x; x = 10 1 (dashed-dotted 
line), x = 1CT 2 (dotted line), x = 10~ 3 (dashed line), x = 10~ 4 (solid line) . 
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- ASYMP 
J k,q 



(39) 



In Fig. [H| we have plotted B,2 )C - Here one observes that this quantity becomes very close to one 
for x < 10~ 2 and Q 2 > 20 GeV 2 which belongs to the region explored by the experiments carried 
out at HERA. This shows that the logarithms mentioned above dominate the structure function 
except in the threshold region given by x ~ 1 and small Q 2 which is characteristic of the EMC 
experiment 0. 

In order to answer the second question whether these logarithmic terms vitiate the pertur- 
bation series one has to resum them in all orders of perturbation theory and show that the 
resummed structure function differs from the one which is computed exactly in fixed order of 
perturbation theory. This resummation procedure is provided by the variable flavor number 
scheme (VFNS) 0. An example of a resummed structure function has been derived in [[L7|]. In 
the case of charm quark production one obtains 



? VFNS 
fc,c 



E 



,,d,s,c 



(/i(4,M 8 )+M4,A* 8 )) 



C2r + /,(4,A* 8 ) 



C k,g 



(40) 
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with the conditions 



'fCI^)^^!, (4i) 




m 2^o "" \m 2 m 2 



lim F™(s, Q 2 , m 2 ) = F fc E c XACT (x, Q 2 , m 2 ) , s = i— E Q 2 . (42) 

s— >4m 2 3? 

The consequence of condition ([y]) is that in the asymptotic regime Q 2 ^> m 2 , F^™ s turns 
into the structure function represented in the four flavor number scheme which contains the 
contribution of light flavors only including the charm quark. The mass singular logarithms, 
occurring in the heavy quark coefficient functions in the three flavor number scheme (see Eqs. 
flT2|), (|T5|)), are shifted to the parton densities defined in the four flavor number scheme appearing 
in Eq. (|40D. The most conspicuous feature is the appearance of the charm quark density which 
is absent in a three flavor number scheme but shows up in the four flavor number scheme. It is 
given by 

/ c (4, f) + / c -(4, /, 2 ) = / g s (3, /i 2 ) ® A% q (j^j + /,(3, /i 2 ) ® A% g {^-\ . (43) 

The operator matrix elements satisfy renormalization group equations which enable us to resum 
all logarithmic terms of the type In fi 2 /m 2 . In order to get the boundary condition (f42|) one 



needs matching conditions for which one can make various choices (see e.g. [||, [17], [fcHJ)- Two 
of them are proposed in [I7|]) ||21|| (BMSN scheme) and in [2!| (CSN scheme). In the former one 
equates 

C NS (^> 4) = C M (^ ), q = u,d,s. (44) 

Both schemes have been calculated up to next-to-next-to-leading order (NNLO) and are com- 
pared in ||I3| with F^ ACT (NLO) in Eq. fl38|). The results are shown in Figs. |8|, ^ from which 



one infers that there is hardly any difference between the two schemes representing VFNS. This 
shows that one can neglect the power contributions 0(m 2 /Q 2 ) in C^ NS which are absent in C^^ q . 
Also the difference between the two versions of VFNS on one hand and the exact NLO approach 
on the other hand is hardly noticeable except in Fig. [5] where in the vicinity of x — 10~ 3 it 
seems that the data are better described by the BMSN and CSN schemes than by the exact 
NLO result. The main conclusion that one can draw from these figures is that the resummation 
effect is very small which means that the so called large logarithms of the type lnQ 2 /m 2 and 
also In /i 2 /m 2 , when /i 2 ~ Q 2 , do not vitiate the perturbation series. 
Summarizing our results we conclude 

1. The past ten years have shown much progress in the computation of higher order corrections 
to heavy flavor production. In particular the results obtained in electro-production agree 
well with the data obtained by the experiments carried out at HERA. 

2. The asymptotic heavy quark coefficient functions can be calculated using operator product 
expansion techniques. The results obtained for the operator matrix elements can be also 
used for processes where the light cone does not dominate the reaction which e.g. holds 



for e + e — > fx + fi 23 
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3. The heavy quark structure function is dominated by the logarithmic terms In Q 2 /m 2 and 
In fi 2 /m 2 provided x and Q 2 are chosen in such a way that they are outside the threshold 
region of the production process i.e. 

s — (1 - x)Q 2 /x > 4m 2 . 

4. In spite of the fact that the logarithms above dominate the structure function they do not 
bedevil the convergence of the perturbation series so that a resummation is in principle not 
necessary. Therefore one can use fixed order (exact) perturbation theory which is simple 
to apply and to interpret in particular if one studies the differential distributions presented 
in Eqs. (|), ©. 
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